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Kinematic dipole signal originates from the motion of the Sun with respect to the CMB rest 
frame. It is a remarkable signal, two orders of magnitude larger than the observed pattern of CMB 
anisotropies. Therefore the dipole subtraction from the main beam is of primary importance. But 
what happens to the dipole signal entering the main spillover (i.e. the most relevant far antenna 
pattern feature, responsible for the bulk of the straylight effect)? We present here an analytical 
and statistical study of this systematic effect for a spinning space-mission like PLANCK * 
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Dipole Straylight Contamination and Low Multipoles 



1. Analytical Description 

The convolution / of the dipole signal with the main spillover 1 can be written as: 

/ = jdClT lm Y{'\e,(p)B SL (d,(p), (1.1) 

where dQ. = ddsinddcp is the solid angle element, 8 G [0,7i] is the colatitude, q> £ [0, 2n [ is the 
longitude, and where the sum on m over —1,0,1 is understood. T\ m are the coefficients of the 
expansion of the dipole 2 on the spherical harmonics basis Y" l {Q,q>) and Bsl(0,(P) is the beam 
response representing the shape of the main spillover in the (6, (p)-plane. In this notation B$l 
is normalized to the whole beam integrated response, dominated by the contribution of the main 
beam J 4n dQ.B ~ / main heam dQ.B ~ 2xoj; where Oj, = FWHM/V81n2. Notice that / is a function 
that depends on the istantaneous orientation of the main spillover. 

Of course eq. ( \l.l\ > is general and exact. Any specific approximation of B$l (inspired by 
experiments) will introduce a certain degree of uncertainty with respect to the real case. Our current 
approximation for Bsl is the Top Hat [jl|, ||] : 

Bsl(Q,9) = /sz.A(0, G ms -A e< , G ms + A 9> )A((p, <jVs-A 9<J (p m . v + A (p> ), (1.2) 

where f$L = p/Q-ms is a constant (p is the ratio between the power entering the main spillover and 
the power entering the main beam, Q. ms is the solid angle subtented by the main spillover) and 
A(a,b,c) = S(a — b)— S{a — c) with S(x) representing the step function (or Heavyside function) 
that takes the value 1 for x > and the value otherwise. Eq. ( |1.2[ ) is just an asymmetric rectangular 
box, in the (0, <p)-plane, centered around the point (d ms , (p ms ) and with sides of length Ag > + Ag< 
and A (p> +A (p< . Notice that the point (6 ms , (p ms ) is the pointing direction of the main spillover. With 
£ and 8 implicitely defined by Aq > = Ag + 8 and A (p> = A v + e where Ae< = Ag and A (p< = A<p, 



we obtain the final expression for the convolution (1.1): 



I TH IfsL = ^ J 2- sin (8 + 2Ag) sin (5 + 2 d ms ) (2A 9 + e) - 
2 V 471 

( A ° + 1 - C ° S {8 + 2 6ms) ^ {5 2 2 ^ ) Re [ T ^ {iP '" s+E/2) } sin + |) ■ d" 3 ) 

where the label TH stands for Top Hat. 

During the rotation of the main beam (assumed, for simplicity, at 90° from the spin axis) in 
each scan circle, the main spillover, if not lying on the spin axis, draws a cone. We define a the 
angle between the main spillover direction and the spin axis 3 . The total signal received by the 
satellite, is the sum of the two following terms: T(d,(p) = Tmb{@, <p) +Isl(@, <P) where Tmb is the 
signal entering the main beam where the dipole has been subtracted away, whereas 1st is the signal 



1 It is the main feature of the antenna pattern far from the main beam. 

2 We use the symbol T( m in order to make clear that the dimensionality is a temperature (°K). Moreover we work in 
a frame centered on the spacecraft with axes pointing fixed (far away) stars. In this frame the dipole vector is constant 
while the direction of the main spillover varies (it rotates of 2k per year in a nominal scanning strategy). 

3 For simplicity we assume here the main spillover centre located on the plane defined by the spin axis and the beam 
centre direction. 



2 



Dipole Straylight Contamination and Low Multipoles 



due to the dipole entering the main spillover. In the case of a PLANCK nominal scanning strategy 
[g], for the first (upper signs) and the second (lower signs) survey, Isl is given by 



I IH (n/2- acos0,«p±7r/2=Fasin0) for < <p< n, 



I 1H (n/2 — acos0,<p=F7r/2±asin0) for % < <p < 2k. 



(1.4) 



The shift in the definition of the arguments in Isl (during either the first or the second survey) 
comes from the fact that when the main beam rotates from North to South the main spillover is 
shifted of — n/2 plus a small correction proportional to a (due to the non perfect alignment of the 
main spillover with the spin axis), while when the main beam rotates from South to North the main 
spillover is shifted of +7l/2 minus a small correction again proportional to a. Notice that (0, <p) is 
the direction of pointing of the main beam. 



2. General Analytical Results 

We set now Tmb = and take into account only the straylight signal, 

T% = J dQI SL (d,(p)Yl m (d,(p), (2.1) 

because we want to analyze only the map due to the Dipole Straylight Contamination (DSC). 
Some general properties can be found PJ|, once performed an expansion for small a: = 
Tfm + a ^im + ^(of 2 )- m particular for an odd number of surveys to the leading order in a 
only the even multipoles survive [i.e. T^j / for even I (and m)] while to the linear order 
only the I = 1 term turns on. For an even number of surveys DSC vanishes except for I = 1, 
= {T^ t] + T^ ld) )/2 = + aT$ + &{a 2 ). These qualitative features have been tested and 
extended by numerical simulations As an example, we show the analytical expressions derived 
for the dipole 

r 1 ^ i = 2 v / 47r/3c 1 a, (2.2) 
Tlfi = (l/2) v / ^/3(±d i +id 2 ){c 2 + c 3 )a, (2.3) 

and the non vanishing coefficients for the quadrupole 

T 2 ±2 = ~ (4/3) 2 V15/327T (di ± 2id 2 ) (c 2 + c 3 ) , (2.4) 

where c x = ^ijAn f SL A sin(2A) T 10 , c 2 = 4^37871 f SL A , c 3 = A^fiJ%Kf SL sin(2A)/2 , d x = 
sin ARe [Tu] , d 2 = sinA/m [Tu] , with Aq = A^ = A and £ = 5 = 0. 



3. Dipole and Quadrupole Analysis 

Section |2| shows that when the main spillover direction is different from the spin axis (i.e. 
then the Dipole itself is modified. This has an impact on the calibration with amplitude 
related to the value of p. For Aq =A,p = A = 7r/10, a = 7r/18 and p = 1/100, we estimate that 
T? L = -2.1nK,Tff = (7.7 ± L(H)flK. This leads to Cf L = 44.5p:K 2 . 
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Figure 1: Left Panel: y = y{F SL ) with B > (lower curves). Solid line is for Cp y = 1500jJ.K 2 , dotted line is 
for Cf 87 = 1000/iK 2 and dashed line is for C$ K¥ = 500jU# 2 . We report also the case of B < (upper curves) 
which does not imply a quadrupole decreasing but only a quadrupole increasing. Right Panel: Distribution 
of the DSC, i.e. C 2 - C% KY . See also the text. 



The next multipole that is modified by DSC is the Quadrupole. We can write the observed 
quadrupole C2 as C2(F SL ) = C|^ F — ^BF$l + ^AF^, where C|^ F is the intrinsic quadrupole (that 
would be observed if the DSC would be absent), B = Re[T n ]Re [T 2 S 2 KY ] +2Im[T n ]Im [T 2 S 2 KY ], 
A = [Re[T n ]] 2 + 4 [Im[T n ]] 2 and F SL = (V5p/3K) [1 + cos A sin A/A]. Depending on the sign of 
B (or in other words on the sign of the real and the imaginary part of T 22 KY ) we have two different 
kinds of behaviour for the observed quadrupole in terms of F$l- It can only increase for B < 
but it can also decrease for B > 0. In the left panel of Fig. [| the two branches are plotted for the 
function y(F SL ) = C 2 (F SL )/C% K¥ where, as an example, it has been set Re [T^ 2 KY ] ~ ±(C% KY /2) 1 / 2 
and where it has been considered Im [Tn] "C Re [Tu]- In this regime Im [T^ 7 ] does not enter in 
the estimate. 

Since the behaviour of the observed quadrupole depends on the sky realization (T^f 7 ), we 
present a statistical analysis aimed at the computation of the probability of increasing or decreas- 
ing of C%. We have made 50000 Gaussian distributed extractions of T^ 7 such that on average 
C|^ y ~ 1000 \iK 2 and we have computed the distribution for the observed quadrupole C2 and for 
the intrinsic quadrupole. In the right panel of Fig. [I] we show the plot of the distribution for the pure 
DSC (i.e. the difference C2 — C|^ y ). We note that the distribution of DSC has a Gaussian profile 
and we compute its mean (202.3 piK 2 ) and standard deviation (404.3 piK 2 ). We conclude that the 
C2 probability of increasing (decreasing) is 69% (31%). 
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